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ABSTRACT 

We numerically integrated the orbits of 1458 particles in the region of the 
classical Kuiper Belt (41 AU < a < 47 AU) to explore the role of dynamical 
instabilities in sculpting the inclination distribution of the classical Kuiper Belt 
Objects (KBOs). We find that the selective removal of low-inclination objects 
by overlapping secular resonances (z/ 17 and acts to raise the mean inclination 
of the surviving population of particles over 4 billion years of interactions with 
Jupiter, Saturn, Uranus and Neptune, though these long-term dynamical effects 
do not themselves appear to explain the discovery of KBOs with inclinations 
near 30°. Our integrations also imply that after 3 billion years of interaction 
with the massive planets, high inclination KBOs more efficiently supply Neptune- 
encountering objects, the likely progenitors of short-period comets, Centaurs, and 
scattered KBOs. The secular resonances at low inclinations may indirectly cause 
this effect by weeding out objects unprotected by mean motion resonances during 
the first 3 billion years. 

Subject headings: celestial mechanics — comets: general — Kuiper Belt — meth- 
ods: n-body simulations — planetary systems: formation 
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1. Introduction 

Most of the mass of the Kuiper Belt between 30 and 50 AU appears to reside in a region 
outside semimajor axis a = 40 AU called the classical Kuiper Belt (Jewitt et al. 1998). 
Early Kuiper Belt surveys found that Kuiper Belt Objects (KBOs) can have surprisingly 
high orbital inclinations (Jewitt & Luu 1995; Jewitt et al. 1996). Recent large-scale surveys 
(Jewitt et al. 1998; Trujillo et al. 2001) have established that the classical KBOs — not just 
the relatively nearby plutinos — frequently have high orbital inclinations (i > 15°). 

The discovery of these high-inclination objects seems to point to unknown processes 
in the primordial solar system. Several mechanisms for "pumping-up" the ancient KBO 
inclination distribution have been investigated, such as resonant encounters (Malhotra 1995; 
Nagasawa & Ida 2000) and perturbations from passing stars (Ida et al. 2000), scattered 
planets (Thommes et al. 2002) and planetesimals (Petit et al. 1999). Large KBOs that may 
have existed long ago during the epoch of planet formation could also have stirred the KBO 
orbits (Kenyon 2002). 

However, theories of the ancient solar system alone can not explain the current KBO 
orbital distribution. Orbit integrations in the Kuiper Belt (Torbett 1989; Torbett & Smolu- 
chowski 1990; Gladman & Duncan 1990; Levison & Duncan 1993; Holman & Wisdom 1993) 
have shown that the dynamical effects of the massive planets on the Kuiper Belt have prob- 
ably removed most of the original KBOs over the lifetime of the solar system. Objects 
removed recently by these processes probably supply today's population of Centaurs, short- 
period comets, and scattered KBOs (Fernandez 1980; Duncan et al. 1988). One might ask 
whether selective removal of objects by interactions with the planets has altered the inclina- 
tion distribution of the surviving population. 

Furthermore, debiased estimates of the KBO inclination distribution (Brown 2001) sug- 
gest that the classical KBOs divide into two populations: one dynamically warm, with a 
typical inclination of i f=a 17°, and one dynamically cold, with a typical inclination % m 2°. We 
wonder whether long-term dynamical interactions helped create the apparent two-component 
distribution. 

To understand the role of dynamical stability in shaping the classical KBO inclination 
distribution and creating Centaurs, comets, and scattered KBOs, we performed a new large- 
scale simulation of the dynamics of KBOs, integrating the orbits of 1458 test particles in 
the region of the classical Kuiper Belt under the influence of gravity from Jupiter, Saturn, 
Uranus, and Neptune. Duncan et al. (1995) performed the most recent previous large-scale 
integration of orbits in the Kuiper Belt. They used an initial grid with fine resolution in 
semimajor axis, but only sampled a few initial inclinations, and did not integrate the orbits 
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of high-inclination objects % > 1° for longer than one billion years. Our work complements 
theirs; at the cost of coverage in semimajor axis, we focus on exploring the consequences of 
high inclination and eccentricity over a four-billion-year period. 

2. The Observed Classical KBO Orbits 

Figure 1 shows the orbital eccentricities and inclinations with respect to the ecliptic 
plane of the 410 known classical KBOs with 40.5 < a < 47.5 AU sorted into seven panels by 
semimajor axis, a, rounded to the nearest AU. We obtained these data from the Minor Planet 
Center's website. 2 Filled symbols indicate orbits for KBOs observed at multiple oppositions. 
The single-opposition orbits, shown as empty squares, generally have accurate inclinations, 
but poorly measured eccentricities; often the Minor Planet Center assumes zero eccentricity 
for single-opposition orbits. 

These classical KBOs have a mean inclination of 6.6°. However, the data reflect a strong 
observational bias against detecting high-inclination KBOs; high inclination objects spend 
less time near the ecliptic plane than low-inclination objects. The figure also reflects severe 
observational biases against detecting distant KBOs. The brightness of a KBO in reflected 
light decreases roughly as 1/r 4 , where r is the KBO's heliocentric distance. This distance 
bias propagates into correlated biases in semimajor axis, eccentricity, and mean anomaly, 
since most KBOs are discovered near perihelion. See Trujillo & Brown (2001) and Gladman 
et al. (2001) for recent discussions of the radial distribution of KBOs. 

Plotted on Figure 1 is a family of curves corresponding to constant values of H = 
a/ GM Q a(l — e 2 ) cosi, the component of the orbital angular momentum (per particle mass) 
in the direction perpendicular to the ecliptic plane. H is largest in the lower left-hand corner 
of each plot (i — e — 0). In a system where all the planets occupy circular orbits in the 
same plane, the orbit of a test particle conserves the component of its angular momentum 
perpendicular to that plane (Kozai 1962; Thomas & Morbidelli 1996). We would expect the 
test particles in our simulation to conserve H and only evolve parallel to these curves in the 
absence of semimajor axis variations if all the perturbing planets had % = e = with respect 
to the ecliptic plane. 

Figure 2 shows the eccentricities and semimajor axes of the known KBOs observed at 
more than one opposition. Two curves on this figure show the eccentricities corresponding 
to perihelion distances of q = 30.1 AU (the Neptune crossing line) and q = 39 AU. Also, each 



2 The Minor Planet Center's website is available at http://cfa-www.harvard.edu/iau/mpc.html. 
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Fig. 1. — Inclinations and eccentricities of the known classical KBOs, sorted by semimajor 
axis rounded to the nearest AU. The vertical lines indicate perihelion distances of 35 and 39 
AU, and the curves show constant values of H = ^jGM Q a(l — e 2 ) cosi, the component of 
angular mementum perpendicular to the ecliptic plane. 
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panel in Figure 1 has vertical grey lines indicating q = 35 AU and q = 39 AU. The dearth 
of objects with q much greater than 39 AU in Figure 2 shows the bias towards detecting 
objects with small perihelia. 

Forty-seven of the classical KBOs with observations at multiple oppositions have orbits 
with perihelia less than 39 AU. These 47 objects have a mean inclination of 13°, much higher 
than the mean inclination of the classical KBOs as a whole (6.6°). Only 11 of these 47 have 
orbits with perihelion q < 35 AU. Duncan et al. (1995) found that particles with q < 35 AU 
generally did not survive one billion years of interaction with the giant planets. 

Grey curves in Figure 2 indicate the boundaries of the 8:5, 5:3, 7:4 and 2:1 mean 
motion resonances (MMRs) with Neptune, using the planar circular restricted three-body 
problem as a model for the motion of small bodies in the field of the Sun and Neptune. We 
numerically computed these boundaries following the technique developed by Ferraz-Mello & 
Sato (1989) to avoid relying on analytic expansions of the disturbing function which invoke 
low eccentricities. Morbidelli et al. (1995) first applied this approach to MMRs in the Kuiper 
belt. The resonance widths calculated this way contain a chaotic zone along the boundary 
at high eccentricities, while the widths calculated by Malhotra (1996) exclude this region. 
This chaos may affect the long-term dynamics and stability of KBOs along the resonance 
boundaries. 



3. Orbit Integrations 

We integrated the orbits of 1458 non-interacting particles for 4 billion years under the 
gravitational influences of Jupiter, Saturn, Uranus, Neptune and the Sun using the MVS 
symplectic integrator from the SWIFT package by Levison & Duncan (1994) based on an 
algorithm by Wisdom & Holman (1991). The initial arguments of pericenter, mean longi- 
tudes, and ascending nodes were chosen randomly from a uniform distribution on the interval 
[0,27r). We used orbital elements referred to the ecliptic plane, rather than the invariable 
plane, so they can easily be compared with observations. We used orbital data for the plan- 
ets from Cohen et al. (1973) and time steps of one year. The integration required donated 
time on 20 different Sun workstations at the Caltech Astronomy Department and at the 
Harvard-Smithsonian Center for Astrophysics. To minimize network traffic, we recorded the 
orbital elements of each particle only once every 10 million years. 

To check that our choice of a one year step size was small enough, we integrated the 
orbits of 100 objects for a billion years in the vicinity of a subtle high-inclination feature 
(a=42, e=0, i=[60,60.02,60.04,. . . ,61.98]; see Section 4.2), using a step size of 1 year. Then 
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Fig. 2. — Eccentricity vs. semimajor axis for the orbits of the known classical KBOs observed 
at more than one opposition. High inclination orbits (i > 15°) are squares; low-inclination 
(i < 15°) orbits are filled diamonds. Grey curves indicate numerically computed boundaries 
of mean motion resonances (MMRs) with Neptune, and black curves show where the peri- 
helion distance, q, is 30.1 AU, i.e., the KBOs become Neptune crossing, (thick curve) and 
where q = 39 AU (thin curve). 
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we repeated the integration using the same initial conditions and a step size of 1/4 of a year. 
Naturally, the final orbits differed due to the non-linearity of the dynamics and the long 
integration times. But the final distributions of the orbits in eccentricity and inclination 
space were indistinguishable. 

Our initial population of objects had semimajor axes of [41, 42, 43,. . . , 47] AU, eccen- 
tricities of [0, 1/30, 2/30,..., 8/30] and inclinations of [0, 1/30, 2/30,..., 17/30] radians 
([0.0000, 1.9099, 3.8197,. . . , 32.4676] degrees). We refer to this grid of 1134 particles, which 
evenly fills a rectangular prism in (a, e, i) space, as our first batch. We also performed some 
exploratory integrations with higher initial inclinations; at 43 and 45 AU, we started par- 
ticles with inclinations up to 1.1666 radians (66.8451 degrees). We removed particles from 
the integration when they encountered Neptune's Hill sphere. Six hundred and twenty-six 
particles from the first batch survived the entire integration. 



3.1. Semimajor Axis Evolution 

Though the particles began with integral values of semimajor axis, after 10 million years 
of integration only a shadow of the discrete initial conditions remained. If we combine all 
the data in the first batch, we find that by t = 10 million years, the particles have spread 
out in semimajor axis into a distribution with half-width at half maximum of 0.4 AU. At 
t = 4 billion years, the core of the distribution has not changed, but the particles in the 
wings are gone. Only 0.1 % of the particles that survived for 4 billion years were transported 
farther than 0.7 AU from their initial semimajor axes. We infer that the observed classical 
KBOs have not strayed far in semimajor axis since the planets settled into their current 
configuration — unless another effect besides long term interactions with the planets disturbed 
their orbits. 

Figure 3 shows the eccentricities and semimajor axes of the particles that survived for 4 
billion years. The eccentricities and semimajor axes have been averaged over the period 3-4 
billion years (100 data dumps). Solid diamonds show the averaged orbits of low inclination 
test particles {i < 15°); empty squares show the averaged orbits of high-inclination particles 
(15° < % < 30°). In the lower right, particles at low e and high a remain clustered, reflecting 
the initial conditions. This zone appears to be stable for the lifetime of the solar system, 
however few KBOs are known with high a and low e orbits, since the perihelia of these 
orbits are far from the sun. To the left of the figure, at a < 44 AU, an inclination-dependent 
instability appears to have removed many of the low-inclination particles, particularly at 
perihelia q < 39 AU. 
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Fig. 3. — Mean eccentricity vs. mean semimajor axis for the particles in the first batch of our 
integration that survived for 4 billion years, averaged over the period t = 3 billion to t — 4 
billion years. Vertical grey lines indicate the nominal locations of mean motion resonances 
(MMRs) with Neptune, and black curves show where q — 30.1 AU (thick curve) and where 
q = 39 AU (thin curve). 
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The grey vertical lines indicate the nominal positions of the 8:5, 5:3, 7:4 and 2:1 MMRs. 
Since the plotted orbits are averaged over the resonant librations, the libration widths of the 
resonances do not appear. Figures 2 and 3 show that objects with small perihelion distances, 
q, often inhabit MMRs, which can protect them from close encounters with Neptune. The 
MMRs at a < 44 AU appear to serve double protective duty; they protect low q objects 
from close encounters, and protect low inclination objects from ejection by the instability 
that removes low inclination particles (see Sections 4.1 and 6). 

3.2. Eccentricity and Inclination Evolution 

Figure 4 shows the orbits of the test particles after ten million years of integration. The 
particles are sorted into seven panels by their initial semimajor axes. Diamonds indicate 
the orbital eccentricity and inclination of the particles that survived the integration for ten 
million years. Thin lines connect the diamonds to the initial orbits. Crosses indicate the 
initial positions of particles that suffered close encounters with Neptune. 

Even after this relatively short period, many of the particles at low inclinations % < 15° 
have been substantially perturbed, their inclinations raised. Particles with inclinations % > 
20° have generally evolved along lines of H =constant, shown in Figure 1. But so far, only 
a few low-inclination particles have been removed. 

Figure 5 shows the state of affairs after four billion years, summarizing the inclination 
and eccentricity evolution of the particles in the same manner as Figure 4. This figure should 
be compared to the current distribution of KBO orbits (Figure 1). By this time, most of the 
low inclination particles whose orbits were greatly perturbed (the diamonds with the long 
lines in Figure 4) have now been ejected. By four billion years, nearly all of the objects with 
initial perihelion q < 35 AU have been removed. This cutoff is not even: several objects 
have been transported into the region q < 35 AU, and at low inclinations % < 10°, the border 
of the unstable region is closer to q = 38 AU (see Figure 1). However, Figure 5 shows a 
systematic removal of high-eccentricity objects, which vividly demonstrates that creating the 
high-inclination KBOs does not require a special mechanism that only increases inclinations 
and not eccentricities. 

When we discuss the relevance of our simulations to the KBO inclination distribution, 
we find it convenient to define two zones in the classical Kuiper Belt, one at 40.5 AU < a < 44 
AU, the other at 44 AU < a < 47.5 AU. This division cuts the distribution of known KBOs 
roughly in half; of the classical KBOs plotted in Figure 1, 213 fall in the first group (84 
with multiple opposition orbits) and 197 fall the second group (92 with multiple opposition 
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Fig. 4. — After 10 million years, the orbits of our test particles have evolved to the locations 
indicated by the diamonds. Short lines connect each diamond to the initial orbit. Crosses 
show the initial locations of particles that have been removed. The particles are sorted into 
seven panels by their initial semimajor axes. 
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Fig. 5. — After four billion years, most of the particles whose orbits dramatically evolved 
in the first ten million years (the diamonds with the long tails in Figure 4) are gone. The 
u 8 , z/ 18 , and Vyj secular resonances have removed some low- inclination particles at initial 
semimajor axes a = 41 and 42 AU. Evolution at high inclinations is mostly parallel to the 
curves of H =constant shown in Figure 1. 
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orbits). The debiased estimate of the radial distribution of KBOs by Trujillo & Brown (2001) 
is also roughly symmetric about 44 AU. As Figure 3 shows, most of the inclination dependent 
structure after 4 billion years is located in the first zone, at a < 44 AU. 

Figure 6 summarizes the inclination distributions of the two zones, using the first batch 
of particles (i < 33°) after 10 million, 1 billion, and 4 billion years. The semimajor axis 
limits are applied when the distributions are evaluated. The 1 billion and 4 billion year 
distributions were averaged over 100 million years (10 data dumps) to cover periodic secular 
motions and provide better statistics. 

After 10 million years (the dashed line), only 9% of the particles have been ejected, 
but many of the low inclination objects have moved to higher inclination orbits. After 4 
billion years, almost half of the particles (46%) are gone, and the region 40.5 AU < a < 44 
AU has been hit the hardest, especially at low inclinations. An approximation to the final 
distribution at 4 billion years (the solid curves in Figure 6) is 

r 0.20 + 0.0017(i - 4°) 2 40.5 AU < a < 44 AU, i < 17.5° 
n(i) « < 0.51 40.5 AU < a < 44 AU, i > 17.5° (1) 

{ 0.7 44 AU < a < 47.5 AU 

4. The Causes of the Instabilities 

As far as we know, all dynamical instability is caused by the overlap of resonances. We 
performed some additional integrations to help understand the dynamics of the unstable 
regions in the classical Kuiper Belt in this context. 

4.1. Secular Resonances at Low Inclinations 

The only notable secular resonances from 41-47 AU are the perihelion and node reso- 
nances with Neptune (u 8 and z/ 18 ) and the node resonance with Uranus (^17), all of which are 
nominally located at a < 41.6 AU (Knezevic et al. 1991). Figure 3 of Duncan et al. (1995) 
shows that these resonances deplete objects at low inclinations at a =41-42 AU. However, 
Figure 4 of Holman & Wisdom (1993) shows that these resonances affect the inclinations 
of particles as distant as a = 50 AU. Our Figure 4 shows in three dimensions a region of 
stirring in inclination and eccentricity that penetrates all of the 7 panels at low inclinations 
and high eccentricities. This region appears to be the broad wing of the zone of influence 
of the overlapped secular resonances centered at 41-42 AU, shown in Figure 4 of Holman & 
Wisdom (1993). 
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Fig. 6. — Histograms of the orbital inclinations of the surviving particles from our first 
batch of integrations (i < 33°) after years (dotted line), 10 million years (dashed line), 1 
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To test this identification, we integrated the orbits of 72 particles with initial a = 
43 AU, e = 0.0, Q = [0°, 54°, 126°, 198°, 270°, 342°], i = [0, 1/50, 2/50, 11/50] radians 
(0°-12.6051°). Ten of these particles were removed before 4 billion years, all from initial 
inclinations < 6/50 radians (6.875°). We plotted pseudo-surfaces of section for all of these 
particles in various coordinate systems to look for librational motion indicating the presence 
of a resonance. 

All of the particles in this group that were ejected shared a common behavior. Figure 7 
shows pseudo-surfaces of section for a typical ejected particle in two coordinate systems: 
(e cos (uj — lvn), e sin (uj — ujn)) and (i cos (Q — Qn), i sin (fl — Qn)), where uj refers to the 
argument of pericenter, Q refers to the argument of ascending node, and the subscript N 
refers to Neptune. The angle uj — ujjy circulates, while the orbit undergoes complex motion 
in the inclination and Q — Q N coordinates, alternating between libration and circulation, 
while the inclination slowly grows. Finding this common behavior indicates that competing 
nodal resonances (vyj and uis) dominate the long-term dynamics of these particles. These 
competing resonances appear to cause the broad swath of destabilization at low inclinations 
and high eccentricities apparent in Figures 3, 4, 5, and 6. 



4.2. The Kozai Resonance 

The Kozai resonance is the tendency of high inclination orbits in the three-body problem 
to undergo secular libration of the argument of pericenter, uj, coupled with oscillations of % 
and e, which conserve H. Recall that we measure an orbit's argument of perihelion from the 
ascending node, so in the solar system, this motion corresponds to an oscillation of uj about 
the invariable plane. Kozai (1962) made the approximation that the planets all orbit in the 
same plane, and calculated that in the solar system, these oscillations set in at inclinations 
of % > 32° for a = 3 AU. Thomas & Morbidelli (1996) calculated that under the same 
approximation, no region of the Kuiper Belt should be subject to the Kozai resonance. 

But Figure 5 shows that the surviving particles around % — 61° show more eccentricity 
evolution than other high inclination particles. The strong inclination dependence of this 
effect and the way the motion appears to conserve H suggest that it is a Kozai phenomenon. 
We integrated the orbits of 24 particles with initial a = 43 AU, e = Vt = at initial uj = 
[0,54°, 126°, 198°, 270°,342°] and initial % = [60°, 61°, 62°, 63°] and examined the evolution 
of their orbits to look for librations, using the coordinates h = ecosw and k = esincj. 
Figure 8 shows three representative pseudo-surfaces of section that illustrate the three types 
of behavior we found in this region. Figures 8a and 8b show the evolution of particles that 
survived the integration, while Figure 8c shows the evolution of one of the two particles 
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(initial u — 0,i — 60° and u = 0,i = 61°) that did not survive. The points in Figure 8a 
fill a roughly circular region, indicating chaotic motion confined to low eccentricities (Henon 
1983). In Figures 8b and c, the motion appears to be a superposition of this chaotic motion 
with libration about uo = 0° or uo = 180°, indicated by the loops. The librational motion 
indicates the influence of a Kozai resonance. 

Occasionally, a particle at low inclinations in a MMR can also undergo Kozai libration 
(Morbidelli 1997). When we plotted pseudo-surfaces of section for all of the 1458 particles in 
the big integration, we found one particle, with initial a = 42 AU, % = 15.278°, e = 0.266666, 
that went into libration about u = 90°. The orbit of this particle, circled in Figures 3 and 5, 
is inside Neptune's 5:3 MMR. The Kozai libration may help the MMR protect the particle 
from doom at the hands of the overlapping secular resonances in this region of dynamical 
space, since the particle survives for the entire 4 billion years of the integration. 

5. The Observed KBO Inclination Distribution 

As we mentioned in Section 2, high-inclination KBOs spend less time near the ecliptic 
plane than low inclination KBOs, so they evade ecliptic plane surveys. Brown (2001) real- 
ized that this observational bias could easily be corrected for the subset of low-eccentricity 
KBOs discovered near the ecliptic plane; for this subset, high inclination objects will be 
underrepresented by the simple factor of sini, assuming no correlation between longitude 
and inclination. The Kozai phenomenon could produce such a correlation, but while Kozai 
libration may be common among plutinos, our integrations in the classical Kuiper Belt pro- 
duced only one Kozai librator at inclinations < 50°. We used the shu debiasing method to 
calculate separate debiased inclination distributions for the inner part of the classical Kuiper 
Belt (40.5 < a < 44 AU) and the outer part of the classical Kuiper Belt (44 < a < 47.5 AU) 
considering only objects discovered at ecliptic latitude < 1°, 209 objects in total. Figure 9 
shows these two distributions, with error bars indicating the Poisson noise due to the limited 
number of objects per bin. 

As Figure 6 shows, 4 billion years of interaction with the massive planets affects the 
inclination distribution of particles in the outer half of the classical Kuiper Belt little, but 
it could substantially alter the inclination distribution of particles with a < 44 AU. Perhaps 
the inclination distribution of the inner half of the classical Kuiper Belt (Figure 9a) once re- 
sembled the inclination distribution of the outer half of the classical Kuiper Belt (Figure 9b). 
We have too few objects to make a detailed comparison, but the distribution in Figure 9b 
clearly has relatively more low-inclination objects than the distribution in Figure 9a. 
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Fig. 7. — Pseudo-surfaces of section for a typical ejected low-inclination particle with initial 
semimajor axis a = 43 AU, eccentricity e = 0. a) A polar plot of eccentricity and uj — uj n b) 
A polar plot of inclination and Q — Qjy- 
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Fig. 8. — Pseudo-surfaces of section for three particles with initial inclination i = 61°, 
semimajor axis a = 43 AU, and eccentricity e = illustrating the three kind of behaviors 
we found near % = 61°. a) Stable chaotic motion confined to low eccentricities, b) Stable 
chaotic motion combined with slow libration around uj = 90° or uj = 270°. c) A particle 
which is ejected under the influence of the Kozai resonance. Note the three different scales. 
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Figure 9a shows a modest peak at % m 8° which piques our curiosity. After only 10 
million years of interaction with the massive planets, we find that a substantial fraction 
of test particles at low initial inclinations % < 5° have been promoted to orbits with higher 
inclinations % ~ 10° (see Figure 4 and Figure 6). This promotion generally fortells an eventual 
close encounter with Neptune. Perhaps the peak at % ~ 8° in Figure 9a represents KBOs 
recently delivered to % ~ 0° by some self-interaction of the KBO disk and promoted to % ~ 8° 
on their way toward being hurled into Neptune's Hill sphere by the overlapped u 17 and vi 8 
secular resonances. 



6. Comet and Scattered KBO Formation 

The classical Kuiper Belt may be the primary source of short-period comets, Centaurs, 
and scattered Kuiper Belt Objects. Levison & Duncan (1997) found that objects from the 
Kuiper Belt that wander into the solar system after a close encounter with Neptune have a 
median dynamical lifetime of only 4.5 x 10 7 years, so perhaps a substantial fraction of short 
period comets and scattered KBOs were classical KBOs less than a billion years ago. To 
explore the origin of short-period comets, Centaurs, and scattered KBOs, we examined the 
particles in our integration which were ejected last from the classical Kuiper Belt. 

Figure 10 shows the initial semimajor axes, eccentricities and inclinations of the particles 
in our integrations that encountered Neptune between t — 2 billion and t — 4 billion years. 
Most (85%) of this population of ejected objects originated on orbits with perihelion distances 
q < 39 AU. Traingles show the initial orbits of particles ejected between t = 2 billion and 
t = 3 billion years. Asterisks show the initial orbits of particles ejected between t = 3 billion 
and t — 4 billion years. 

We were surprised to find that the objects ejected from our simulation between t = 3 
billion years and t — 4 billion years came exclusively from high initial inclinations (i > 15°). 
Objects in our first batch encountered Neptune at a rate roughly proportional to t~ l for the 
first 1 billion years, in agreement with Holman & Wisdom (1993) and Duncan et al. (1995). 
However, during the last billion years of the simulation, not one particle from low inclinations 
(i < 15°) encountered Neptune, while 14 particles from higher inclinations (15° < i < 33°) 
encountered Neptune during this interval. We would expect ~ 10 encounters for the low 
inclination group during that period if the rate remained proportional to t~ x . 

Naturally, the overlapping secular resonances that create the large region of instability 
at low inclinations throughout the classical Kuiper Belt affect the rate of Neptune encounters 
for low-inclination objects. These secular resonances remove low-inclination objects early on, 
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Fig. 10. — Initial orbits of particles that are removed late in the integrations when they enter 
Neptune's Hill sphere. Most of these particles originate on orbits with perihelia between 35 
and 39 AU (the vertical grey lines). Particles removed between t = 2 billion and t = 3 
billion years (the triangles) come from all inclinations. Particles removed between t = 3 
billion years and t = 4 billion years originated only at high inclinations (i > 15°). 
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so relatively few remain at t > 3 billion years to encounter Neptune; in the region q < 39 
AU almost twice as many high inclination (15° < % < 33°) objects survive for 3 billion 
years as low inclination objects (i < 15°). However, this selective removal of low inclination 
objects does not explain why not a single low inclination object in our simulation encountered 
Neptune from t = 3 billion to t — 4 billion years. Enough low inclination objects remain at 
t — 4 billion years to provide a significant flux of Neptune encounters. 

But the secular resonances do not act alone. MMRs can protect the low-inclination 
objects from destruction at the hands of the overlapping secular resonances. The secular 
resonances may weed out all but the MMR-protected low-inclination objects in the first 3 
billion years, so after that time, the only available unprotected objects for them to eject have 
high inclinations. 

Figure 3 shows that at a < 44 AU, all but a few particles with q < 39 AU that survive 
at low inclinations for 4 billion years (Figure 3) inhabit MMRs. High inclination particles 
in this zone are spread more evenly; some high inclination particles lie in MMRs, but many 
solid squares fall between the MMRs. Figure 2 shows that the real KBOs mimic the behavior 
of the test particles which survive 4 billion years; 12 out of the 15 low inclination (i < 15°) 
KBOs at a < 44 AU, q < 39 AU with multiple-opposition orbits have observed semimajor 
axes and eccentricities within or near the numerically calculated boundaries of MMRs. The 
observed high-inclination KBOs are spread more evenly in semimajor axis. 

Our simulation suggests that todays short-period comets, Centaurs, and scattered KBOs 
originate in the high-inclination KBO population. Therefore, short-period comets, Centaurs, 
and scattered KBOs should have a color distribution like the high-inclination KBOs when 
they leave the Kuiper belt. Trujillo & Brown (2002) found that low inclination KBOs are 
uniformly red (B — R « 1.5-2.2 for i < 12°), while high-inclination KBOs range from red to 
blue (B-Rm 1.0-2.0 for i > 12°). 

Jewitt (2002) compared the colors of KBOs, Centaurs, comet nuclei, and candidate dead 
comets, and found that the Centaurs in his study had a wide range of colors. This range 
appears to be consistent with the color range of the high-inclination KBOs in (Trujillo & 
Brown 2002), in agreement with our dynamical result. More color measurements of Centaurs 
may determine whether the color distribution of centaurs matches the color distribution of 
the high-inclination KBOs in detail. Jewitt 's samples of comet nuclei and candidate dead 
comets had no objects as red as the red high-inclination KBOs; our study offers no solution 
to the problem of this missing ultrared matter. 
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7. Conclusions 

Dynamical erosion is not the only long-term effect in the classical Kuiper Belt; collisions 
may shape the distributions of KBO orbital elements (Stern & Colwell 1997; Kenyon & Luu 
1998, 1999; Durda & Stern 2000). But as the collision rate and the mass of the Kuiper 
Belt decreased with time, the dynamical effects shown by our integrations must have begun 
to dominate the shaping of the classical Kuiper Belt. We found that interactions with the 
massive planets preferentially deplete low-inclination objects in the inner half of the classical 
Kuiper Belt, raising the mean inclination of the population of surviving objects. However, 
the outer half of the classical Kuiper Belt (a > 44 AU, q < 39AU) could have retained its 
ancient inclination distribution. 

We also found that objects ejected from the classical Kuiper Belt during the last 1 
billion years of our integration primarily come from initial orbits at high inclinations and 
eccentricities. This effect may be caused indirectly by vyj and z/ 18 , which remove many low- 
inclination particles at early times and weed out those unprotected by MMRs. This finding 
implies that Centaurs, short-period comets and scattered KBOs should have an initial color 
distribution like that of the high-inclination KBOs. Measuring of the colors of more Centaurs 
could test this hypothesis. 

Perhaps the high-inclination KBOs separated from the primordial KBOs objects during 
a special event early in the lifetime of the solar system. If today's active comets arise 
mainly from high-inclination KBOs, we have this event to thank for their inner solar system 
company. 
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